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We propose and analyze theoretically an experimental setup for detecting the elusive Majorana 
particle in semiconductor-superconductor heterostructures. The experimental system consists of 
one-dimensional semiconductor wire with strong spin-orbit Rashba interaction embedded into a su- 
perconducting quantum interference device. We show that the energy spectra of the Andreev bound 
states at the junction are qualitatively different in topologically trivial (i.e. not containing any Majo- 
rana) and nontrivial phases having an even and odd number of crossings at zero energy, respectively. 
The measurement of the supercurrent through the junction allows one to discern topologically dis- 
tinct phases and observe a topological phase transition by changing the in-plane magnetic field or 
the gate voltage. The observation of this phase transition will be a direct demonstration of the 
existence of Majorana particles. 
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The Majorana fermions were envisioned by Ettore Ma- 
jorana [1 in 1937 as fundamental constituents of nature. 
Majorana particles are intriguing and exotic because each 
Majorana particle is its own antiparticle unlike Dirac 
fermions where electrons and positrons (or holes) are dis- 
tinct. Recently, the search for Majorana fermions has fo- 
cused on solid state systems where many-body ground 
states may have fundamental quasiparticle excitations 
which are Majorana fermions [2 . Although the emer- 
gence of Majorana excitations, which are effectively frac- 
tionalized objects (anyons) obeying non-Abelian anyonic 
statistics rather than Fermi or Bose statistics [3 , in solid 
state systems is by itself an extraordinary phenomenon, 
what has attracted a great deal of attention is the pos- 
sibility of carrying out fault tolerant topological quan- 
tum computation in 2D systems using these Majorana 
particles [4 . Such topological quantum computation, in 
contrast to ordinary quantum computation, would not 
require any quantum error correction since the Majo- 
rana excitations are immune to local noise by virtue of 
their nonlocal 'topological' (TP) nature [3l[4]. The direct 
experimental observation of Majorana particles in solid 
state systems would therefore be a true breakthrough 
both from the perspective of fundamental physics of frac- 
tional statistics in nature and the technological perspec- 
tive of building a working quantum computer. It is there- 
fore not surprising that there have been several recent 
proposals for the experimental realization of Majorana 
fermions (MFs) in solid state systems 0-^. 

In this Letter, we propose and validate theoretically a 
specific experimental setup for the direct observation of 
MFs in one of the simplest proposed solid state systems - 
ID semiconductor/superconductor heterostructure based 
quantum wires. This particular heterostructure consist- 
ing of an ordinary superconductor {e.g. Nb) and a semi- 
conductor with strong spin-orbit coupling {e.g. InAs) as 



proposed originally by Sau et al. [F and expanded by 
Alicea [7 , is simple and does not require any specialized 
materials for producing Majorana modes. The supercon- 
ductor (SC) induces superconductivity in the semicon- 
ductor (SM) where the presence of spin-orbit coupling 
leads to the existence of MFs at the ends of the wire. 
We show that in a suitable geometry (see Figjl]) the SC 
state in the semiconductor undergoes a phase transition, 
as the chemical potential or magnetic field is tuned, from 
a superconducting state containing Majorana modes at 
the junction to an ordinary SC state with no Majorana 
modes at the junction. We establish that such a tran- 
sition is indeed feasible to observe in the laboratory in 
semiconductor nanowires, showing in the process how 
one can experimentally discover the Majorana mode in 
the SM/SC heterostructure. 

Specifically, we consider here ID InAs nanowire 
proximity-coupled with an s-wave superconductor {e.g., 
Nb or Al). InAs nanowires in proximity to Nb and Al 
have been studied experimentally [9 and are known to 
form highly transparent interfaces for electrons allowing 
one to induce a large SC gap Aq in InAs (Aq ^ ANb^ 
15K) [10 . Moreover, in this quasi ID geometry (see 
Fig{T]3) the in-plane magnetic field can open up a gap 
in the spectrum at zero momentum and eliminate fermion 
doubling. Because of the vast difference in the g-factor 
for Nb ^Nb^^l and InAs ^inAs^35 [11 , the in-plane mag- 
netic field Bx ^O.IT can open a sizable Zeeman gap in 
InAs {Vx ^ IK) without substantially suppressing SC in 
Nb {H^^^0.2T). The nanowire can be gated [9 allowing 
one to control chemical potential in it. Thus, the current 
proposal involves a simple architecture and yet preserves 
the parameter phase space flexibility, which puts the re- 
alization of MFs in the SM/SC heterostructure within 
the experimental reach. 

We show below that the supercurrent through SM/SC 
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FIG. 1: (Color online) (a) Top view of SM/SC heterostructure 
embedded into small-inductance SC loop, (b) Side view of the 
SM/SC heterostructure. The nanowire can be top-gated to 
control chemical potential. Here we assume L<^^ and Li^^ 
with ^ being the SC coherence length, (c) Proposed read-out 
scheme for the Andreev energy levels. Inductively coupled 
rf-driven tank circuit allows time-resolved measuring of the 
effective state-dependent Josephson inductance [19]. 



heterostructure exhibits unusual behavior due to the 
presence of MFs in the system. In particular, the spec- 
trum of Andreev states has an odd number of crossings 
at E = in the TP phase {Co= IJ^^^Al-V^<0 with /i be- 
ing chemical potential) whereas in the TP trivial phase 
(Co > 0) the number of crossings is even. Odd number 
of crossings is associated with the presence of MFs in 
the system leading to 47r-periodic Andreev energy spec- 
trum [8]. Thus, this difference in the spectrum allows 
distinguishing TP and conventional SCs. The remark- 
able feature of the present proposal is that by changing 
Bx or /i across the phase boundary between TP trivial 
and nontrivial superconducting phases (Co = 0) one can 
contrast different qualitative dependence of the Andreev 
energy spectrum on magnetic flux <I> through the SQUID. 

Theoretical model. We consider an infinite (Li ID 
semiconducting wire embedded into SQUID, see Fig.[l^. 
The Hamiltonian describing the nanowire reads {h = 1) 



poo / q2 

Ho=J dxiljl{x)\-^^-/j.^iaaydx^Vx(J, 



^.'{x), (1) 



where m*, ji and a are the effective mass, chemical po- 
tential and strength of spin-orbit Rashba interaction, re- 
spectively. In-plane magnetic field leads to spin split- 
ting Vx=gsM/J^BBx/2. The radius of the wire R is small 
compared to the Fermi wavelength R<Xf so that there 
is a single ID mode occupied. Because of the proximity 
effect between SM and SC (see Fig.[l]3), Cooper pairs can 
tunnel into the nanowire. These correlations can be de- 
scribed by Hsc = J^oo^x ^A(x)?/^j(x)?/^|(x) + /i.c.^ . Here 

A{x) is the induced pairing potential in the nanowire 
A{x) = Ao6(x-L) + Aoe^^6(-x-L) with cp being the 
phase of the order parameter. 

One can recast the full Hamiltonian H = Hijj-Hsc in the 
dimensionless form by introducing rescaled coordinates 
x = m*ax and energies E = E/m*a'^. The BdG equations 



then become ^BdG^(^) = E'^{x). Using the convention 
for Nambu spinors ^(x) = {u^{x)^ui{x),vi{x)^—v^{x)) 
the BdG Hamiltonian reads 



HbcLG 



1 



dx^i(Tydx-fl Tz^VxCJx 



(2) 



- AQ{x — L)rx -\- AQ{—x — L) (cos (prx-\- sin (pry) . 



The solution of the BdG equations supplemented with 
appropriate boundary conditions yields the Andreev 
spectrum in the junction. It is useful to solve for the 
energy at = tt. At this point the profile of the order 
parameter in the limit of L <C ^ forms a domain wall, 
which under certain conditions can host a pair of Ma- 
jor ana bound states [6 . To demonstrate this we inves- 
tigate the existence of zero-energy solution by solving 
^BdG^o(^) =0. At (p=7r^ BdG Hamiltonian ([2| is real 
and, thus, one can construct real Nambu spinors ^o(^)- 
According to the particle- hole symmetry if ^o(^) is a 
solution, then ayry'^o{x) is also a solution. This im- 
poses the constraint on the spinor degrees of freedom: 
v^/l{x) = Xu^/^{x) with A = ±1. Thus, the 4 x 4 BdG 
Hamiltonian can be reduced to 2 x 2 matrix: 



2^x 



Vx^XA{x)^dx\f u^{x)\_ 



Vx-XA{x)-dx 



u^{x) 



= 0. (3) 



One can seek solutions of Eq. ([3| in the form u^i^{x)(xe^^ 
and require solutions for x ^ to be normalizable. Let 
us concentrate on the x>0 case. Then, the characteristic 
equation for z following from Eq.(|3| reads 

z'^+4(/i+l)z^ + 8AAoz+4Co = with Co = M^+Ag-V;^ 

(4) 

The roots Zi of the above quartic equation with real coef- 
ficients should satisfy the following constraints: nt=i ~ 
4Co and Xl^i = 0- If ^ii I'^^i ^^^d Co > 0, these 

constraints are satisfied only when the number of solu- 
tions with Re[z] ^ is the same. If Eq.Q has at least 
one complex solution zi =a-\-ib^ then there is another so- 
lution Z2 = a—ib. Since the other two solutions are given 
by the quadratic equation, one can express these roots in 
terms of a and b: ^3,4 = —a± y^^^^^^^Cq/I^^^F^- Given 
that |Re[v'a2-4Co/(a2 + 62)]| ^ |^| fo^ Co > 0, there are 
two solutions with Re[z] ^0, respectively. Different val- 
ues of A change the sign of a, and this conclusion is valid 
for both channels A = il. Thus, when Co > there are 
two exponentially decaying solutions for x ^ yielding 
4 coefficients to match. Since the number of constraints 
(4 from boundary conditions and 1 from normalization) 
is larger than the number of linearly independent coeffi- 
cients, there are no zero energy solutions for Co > 0. On 
the other hand, similar analysis for Co < always yields 
three roots with Re[z] <0 either in A = l or A = — 1 chan- 
nels resulting in six coefficients to match. Therefore, in 
this case there is a pair of zero-energy Majorana states. 
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FIG. 2: (Color online) Andreev energy spectrum in SM/SC 
heterostructure for the junction with L ^ 0. (a) Energy 
spectrum in TP trivial (dashed line: \4 = 0.75) and nontrivial 
(solid line: Vx — 1.25) states. The two TP distinct phases 
differ by having even and odd number of crossings, respec- 
tively, (b) Schematic plot of the Josephson current through 
the junction carried by Andreev states: light (red) and dark 
(blue) lines describe Josephson current in TP trivial and non- 
trivial phases, respectively, (c) and (d) The evolution of An- 
dreev energy spectrum with chemical potential, (c) The spec- 
trum in TP nontrivial phase. The dashed (red) line is a fit to 
zbcos((/?/2) function, (d) The spectrum in TP trivial phase. 
There is no crossing at — i^. 



At Co = 0, there is a solution with z = 0, which corre- 
sponds to the closing of the SC bulk excitation gap [6]. 
Therefore, the condition Co = gives the phase boundary 
between TP trivial and nontrivial SC phases [12 . 

Andreev spectrum as a function of magnetic flux ^ 
can be obtained by solving BdG equations defined by 
Eq.Q in the limit of L ^ (describing L <C ^ case) 
and matching the boundary conditions ^(0~) = ^(0+), 
9aj^(0~)=9a^^(0+). The algebra is not particularly en- 
lightening so we present here numerical results shown in 
Fig. [2) which are consistent with above analytical consid- 
erations. The characteristic signature of the TP nontriv- 
ial phase is the presence of odd number of crossings in the 
Andreev spectrum in contrast with the TP trivial phase 
where number of crossings is even as required by 27r- 
periodicity of the BdG Hamiltonian, see Fig. [2^. Indeed, 
in the absence of the degenerate TP sectors, upon the ad- 
vance of the SC phase by 27rthe system returns to the 
same state. It is well-known that in S C- normal- met al-SC 
heterostructure the spectrum of spin-degenerate Andreev 
states is ^((^) = Ao\/l — sin((^/2)^ [14], where D is the 
interface transparency. The presence of weak spin-orbit 
interactions leads to the degeneracy splitting of Andreev 
levels [15]. In the TP trivial phase, which is adiabatically 
connected to T^^O limit, we obtain similar results, see 
Fig. [2]i. In contrast, as shown in Fig. [2]3 Andreev spec- 
trum in the TP nontrivial phase is strikingly different. 
This difference is related to the presence of the Major ana 



zero-energy modes in the system at (/? = 7r. The quantum 
phase transition between these two phases is called topo- 
logical because it occurs without any qualitative changes 
of the local order parameter. The two phases are distin- 
guished by the topological order associated with the pres- 
ence of Majorana zero-energy modes. The TP quantum 
phase transition occurs when A^, which is proportional 
to the quasiparticle bulk gap, becomes zero bringing a 
continuum of gapless states at £^ = 0. This phenomenon 
is generic and applies also to Majorana bound states in 
the vortex cores. The topological reconstruction of the 
fermionic spectrum cannot occur adiabatically and re- 
quires the nullification of the bulk excitation gap [ 13] . 
Looking at Figs. [2^ and |2|3, one can see the evolution 
of the Andreev energy spectrum with the magnetic field: 

= 2, 1.25, 0.75, which supports above arguments. Also, 
Figs.jSJ) and|2|i show the evolution of the spectrum with 
the chemical potential. 

We note that the position of the zero-energy crossing is 
not universal and can be shifted by adding a weak pertur- 
bation, e.g. ByCFy. However, the crossing itself is robust 
and is protected by particle-hole symmetry. Indeed, the 
eigenstates with are related by particle-hole symme- 
try = 0^-£;(x), where 6 = (JyTyK with K being 
the complex conjugation operator. One can show using 
the property cFyTyHcFyTy = —H^ that matrix elements 
(^|i^6|^) = -(^|i^6|^) =0, and thus, the crossing is 
protected against any perturbations as long as the bulk 
gap is preserved. Another elegant way of demonstrat- 
ing the robustness of the crossing point was suggested in 
Refs. [81117]. At E = and (p = 7r one can introduce two 
MF operators 71,2 =7i 2- Then, the low-energy Hamilto- 
nian around cp ^ tt can be written as H = i2e{(p)jij2- By 
introducing the Dirac fermion operators = 71+^72 and 

= 71 — Z72 , one can rewrite the Hamiltonian above as 
H = s{(p){c^ c—^), from which it follows that the states 
and "^-E have different fermion parity. Thus, as long as 
fermion parity is locally conserved the matrix elements 
between the states "^e and "^-e are zero. 

Two (or even number) crossings in the Andreev spec- 
trum as in Fig. |2]i are not generally protected. We 
have studied the robustness of even and odd crossings 
numerically by adding the impurity-scattering potential 
U{x) = UoS{x-L) into Eq.Q. As shown in Fig. ^ 
impurity scattering opens up a gap in the spectrum in- 
dicating that crossings in the TP trivial phase are not 
robust. In contrast, impurity scattering does not affect 
the crossing in the TP nontrivial phase (see Fig. |3^) . We 
also considered finite-size heterostructure L = 3, where 
the spectrum has excited Andreev states. As shown in 
Fig. [3^, the crossing at zero energy is robust while other 
crossings are not. 

The experimental system shown in Fig. can be 
viewed as two Majorana quantum wires [8 coupled by 
tunneling through the junction. Indeed, consider a 
SM wire of length Li at x > 0. One can diagonal- 
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FIG. 3: (a) Andreev spectrum for a finite-size junction L = 3 
in a TP nontrivial phase. Here /i = 0, A = l, 14 = 2 and 
Uo/a = 1. (b) Andreev spectrum in TP trivial phase for 
= 5, A = 1, T4 = 2, L < 1 and Uo/a 1. 

ize the single-particle Hamilt onian Q and find eigenval- 
ues e±{p) =p^/2m*— /i ± and eigenvectors 

^±(p) = 7^(±(^^+^<^P:r)/\A^+^^'l)^- Assuming 
that only lowest band €-{p) is occupied, the full Hamil- 
tonian H can be projected to the lowest band yielding 

p 

where the order parameter A_{p) = iapxAo/ y/ a'^p'^ + 
has p-wave symmetry. Thus, the present problem is iso- 
morphic to the Majorana wire considered by Kitaev [8]. 
As long as Li ^ ^, tunneling amplitude of MFs be- 

tween the ends of the wire vanishes toce ^ and differ- 
ent fermion parity ground states are almost degenerate. 
When two Majorana wires are brought together as shown 
in Fig. [TJd, zero-energy Majorana modes at the junction 
are hybridized yielding the spectrum shown in Fig. 
The presence of MFs in the system can be characterized 
by Z2 topological invariant M{Ho) = (-l)^(o)-^(A) [^^ 
where iy{0) and v{A) are the number of negative eigen- 
values of Hq at p = 0. A, respectively. Here A is the mo- 
mentum at the edge of the Brillouin zone. The difference 
z^(0)— i^(A) counts the number of bands (mod 2) cross- 
ing Fermi level on the interval (0,A). In weak pairing 
limit Ao <^fi,Vx, this definition of TP trivial (M = +1) 
and nontrivial phases (M = — 1) is consistent with exact 
results for this model discussed after Eq. Q. 

The difference in Andreev spectrum should be de- 
tectable by various experimental techniques. In partic- 
ular, the Josephson current in L ^ limit is given by 
= -f [H [IS]. The energy Ei^2{^) close to 

(p = 7r is well approximated by ±cos((/:?/2). Thus, the cur- 
rent carried by the quasiparticle state n at = tt is max- 
imum in the TP nontrivial phase in contrast to the TP 
trivial case where /n = 0? see Fig. (2}). This phenomenon 
was dubbed fractional Josephson effect O [161 (HI- 
reality, however, there are processes changing fermion 
parity, and current will fluctuate between I± = ±7 with 
switching time r. Such processes were studied in the 
context of SC qubits fT8|, where the fermion parity 



switching time r was measured experimentally yielding 
r > 1ms at T = 20mK in Al. The random telegraph 
signal of Josephson current can be measured by induc- 
tively coupling the SQUID to the rf-driven tank circuit 
(see Fig. ^) and monitoring in real time the impedance 
of the circuit, which depends on effective Josephson in- 
ductance Lj^((/p )= ^2- ^-^^ [19^- For typical parameters 
of InAs m* ~ 0.04me, a ~ O.leVA corresponding to the 
leuG^th scale -^^lOOnm and Vx^lK, An^lK, the criti- 

cal current - lOnA and ^L^™\(/p)| -10- lOOnH. The 
Josephson inductance Lj{(p) can be probed by small- 
amplitude phase oscillations with the frequency uj sat- 
isfying cj^; A^~ lOGHz (for adiabatic approximation to 
hold) [20] and uj:^l/r (to resolve current fluctuations). 
Thus, this experimental technique can be used to distin- 
guish the Andreev spectrum in TP distinct phases and 
observe the phase transition we predict. 

This work is supported by DARPA-QuEST and JQI- 
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